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We compute the renormalized expectation value of the square of a quantum scalar field on a 
Reissner-Nordstrom-de Sitter black hole in which the temperatures of the event and cosmological 
horizons are equal ('lukewarm' black hole). Our numerical calculations for a thermal state at the 
same temperature as the two horizons indicate that this renormalized expectation value is regular on 
both the event and cosmological horizons. We are able to show analytically, using an approximation 
for the field modes near the horizons, that this is indeed the case. 
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PACS numbers: 04.62+v,04.70.Dy 

I. INTRODUCTION 

The renormalized stress-energy tensor (RSET) 
(T^ren is an object of fundamental importance in 
quantum field theory in curved space-time, since it 
governs, via the semi-classical Einstein equations 

G>„ = 8TrG(Tfj, v ) Icn , 

the back-reaction of the quantum field on the space-time 
geometry. Computing the RSET on a particular space- 
time background is a complicated process [H, @, d, Q , and, 
for a quantum scalar field (f>, it is informative to study 
first the renormalized vacuum polarization (0 2 ) rcn , which 
is considerably easier to compute and has many of the 
same features as the full RSET. For example, although 
(</> 2 )ren is a scalar object and hence cannot distinguish 
between future and past event horizons, nonetheless, if 
it diverges at a horizon for a particular quantum state, 
then it is likely that the RSET also diverges there. 

The vacuum polarization (</> 2 ) ren has been extensively 
studied by many authors for various black hole back- 
grounds, beginning with computations by Candelas [f| of 
(</> 2 )ren for a massless, conformally coupled, scalar field 
on the event horizon (and at infinity) for a Schwarzschild 
black hole. This was subsequently extended to the whole 
of the Schwarzschild geometry in Refs. [6( (exterior to 
the event horizon) and [7J (interior to the event horizon) . 
The corresponding calculation for massive scalar fields 
was done by Anderson Other examples of calcula- 
tions of (0 2 )ron on more general black hole space-times 
can be found in [?|. As well as exact, numerical calcula- 
tions, a number of approximation schemes have been de- 
veloped for various types of space-times, both for (4> 2 ) len 
[3 and the RSET Q. 

Of the three standard vacua for quantum fields on 
black hole space-times (Hartle-Hawking Unruh [l3[ 
or Boulware [3), the Hartle-Hawking state has received 
the most attention in the literature. This is because the 
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Hartle-Hawking vacuum possesses the most symmetries 
(for example, time reversal symmetry) and so is more 
straightforward to calculate. Expectation values of ob- 
servables in this state are also expected to be regular 
on the event horizon of a black hole. Once the renor- 
malized expectation values in a particular state (say the 
Hartle-Hawking state) have been computed, it is much 
easier to compute renormalized expectation values in an- 
other state because the difference between expectation 
values in two states does not require renormalization. 
A general method for computing (</> 2 ) ren on any static, 
spherically symmetric, black hole space-time for a state 
at either a fixed non-zero temperature or zero tempera- 
ture was developed by Anderson [l5j], and subsequently 
extended to calculations of the RSET [l[ (see also [l6| 
for an application of this type of approach to the space- 
time outside a star). This is the method we shall adopt 
in this paper. This approach works for computations of 
the Hartle-Hawking (non-zero temperature) and Boul- 
ware (zero temperature) states, but not for the Unruh 
vacuum. 

Most of the work on (0 2 ) ron and the RSET to date 
has focussed on asymptotically flat black holes, such as 
Schwarzschild and Reissner-Nordstrdm. However, there 
are many interesting features of quantum field theory on 
asymptotically de Sitter black holes [13], when we have 
both a black hole event horizon and a cosmological hori- 
zon. The Kay-Wald theorem [l8j states that there is 
no thermal state on Schwarzschild-de Sitter black holes 
which preserves all the symmetries of the metric (includ- 
ing time-reversal symmetry) and which is regular on both 
the event and cosmological horizons. Therefore there is 
no equivalent of the Hartle-Hawking state for this geom- 
etry. This may be understood heuristically as follows. A 
thermal state is expected to be regular at a horizon if 
the temperature of the state matches the temperature of 
the horizon. However, for Schwarzschild-de Sitter black 
holes, the temperatures of the event and cosmological 
horizons are not equal and therefore any thermal state 
cannot match both temperatures. If we consider the more 
general case of Reissner-Nordstrom-de Sitter black holes, 
then it is possible for the event and cosmological horizons 
to have the same temperature [l9[ (these are known as 
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'lukewarm' black holes 20]). Therefore a natural ques- 
tion is whether the thermal state, constructed on these 
black holes at the natural temperature, is regular on both 
the event and cosmological horizons. A simple calcula- 
tion for two-dimensional black holes [2l| shows that this 
is the case, but for four-dimensional black holes there 
is an unknown function in the RSET which can only be 
found by direct calculation. We will return to the compu- 
tation of this function in the near future 22], but in the 
present paper we will focus on the simpler computation 
of (0 2 ) rcn for this state. 

The outline of this paper is as follows. In section |TT] 
we review the method of [l[ to calculate (0 2 ) rO n using 
point-splitting. As we are considering a thermal quan- 
tum state, the Euclidean Green's function will be used. 
Using temporal point-splitting, even with the points split 
the Green's function contains apparent divergences. We 
present in section III Bl a new approach to regularizing 
these divergences, based on dimensional reduction. Once 
the point-split Green's function is manifestly finite, the 
renormalization procedure is then relatively straightfor- 
ward, following [lj. The renormalized expectation value 
(0 2 )rcn is computed numerically for lukewarm black holes 
in section Hm Our numerical results indicate that (</> 2 ) ren 
is indeed regular on both the event and cosmological hori- 
zons. This is proved analytically in section HVl in which 
we use an approximation for the field modes near the 
horizons, developed in Ref. [23] for the asymptotically 
flat Reissner-Nordstrom case. Finally, our conclusions 
are presented in section |Vj Throughout this paper the 
(Lorentzian) metric has signature (—,+,+,+) and we 
use units in which 8ttG = h = c = kn = 1. 



II. 



GENERAL METHOD TO CALCULATE <0 2 ) r 



A. Point-split Green's Function 

We use the method of Ref. P, to compute the 
renormalized expectation value (0 2 ) rO n- In this section 
we just briefly outline the key steps in the construction, 
further details of which can be found in Refs. [H.[l5l]. 

We begin with a scalar field <j) with mass m and cou- 
pling £ to the scalar curvature R, satisfying the field equa- 
tion 



hole geometries, including, in particular, the Reissner- 
Nordstrom-de Sitter black holes for which 



f(r) = 1 + — - 

r r z 3 



(2.3) 



where M, Q are related, respectively, to the mass and 
charge of the black hole [24[ and A is the (positive) cos- 
mological constant. Further details of the space-time ge- 
ometry in which we are particularly interested, the luke- 
warm black holes, will be given in section IIII Al 

We are interested in computing (</> 2 )_ren for a thermal 
state at a temperature T. We follow [l| and use a Eu- 
clidean space approach. By defining the Euclidean time 
r as r = it, the metric (12.21) becomes 



>dip 2 . (2.4) 



ds 2 = fi^d^+fir)- 1 dr 2 +r 2 dd 2 +r 2 



We define G E (x; x') to be the Euclidean Green's function, 
which satisfies the equation [l[ 

[Va,„V£ - m 2 - £R] G E (x;x') = -g~* {x)8\x, x'), 

where the covariant derivative is now taken with respect 
to the metric (|2 .4[) . Using the method of point-splitting, 
the (unrenormalized) expectation value (</> 2 ), 
by the following limit: 



is given 



lim G E (x; x) 



(2.6) 



The stress-energy tensor is calculated from derivatives of 
G E {x-x') Q. 

For a thermal state at a temperature T, the Euclidean 
Green's function G E (x;x') is periodic in r — r' with pe- 
riod T . In this case the scalar Euclidean Green's func- 
tion takes the form [l], Il5j 

rp 00 

G E {x;x') = — } exp [iu> (r - r')] 



^(2^+l)P,(cos 7 )T^(r,r'),(2.7) 



where u> = 2imT, and Pi is the usual Legendre function, 
with 



(V^V 1 - TO" 



-£R)(/,= 0. 



(2.1) 



In this paper we consider a black hole space-time with 
metric, in Schwarzschild-likc co-ordinates, given by 



ds 2 = -f(r) dt 2 



firy'dr 2 



r 2 d6 2 



r 2 sin 2 9 dip 2 , 



(2.2) 

where the metric function / depends only on the radial 
co-ordinate r. The metric (|2.2j) is not the most gen- 
eral spherically symmetric black hole metric, and, in- 
deed, the method of Ref. [l[ is developed for more gen- 
eral spherically symmetric metrics. However, the metric 
is sufficiently general to cover a broad class of black 



cos 7 = cos 6 cos 6' + sin 9 sin 9' cos (<p — cp ) , 
and T^i satisfies the differential equation 

'2/ , df\ dT u 



1 d 2 T 
—S(r-r') = f— — 
r 2 ( ' J dr 2 



T 



r 

e(£ + i) 



dr ) dr 
- 2 1 ZR 



m 



(2.8) 



T 



(2.9) 



The differential equation (|2.9p arises from separating the 
wave equation (|2.5[) on the background metric 
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We define functions p^e and q u e, as solutions of the 
corresponding homogeneous differential equation 



d 2 T 
dr 2 

f 



'2/ + df \ dr 

r dr 
1(1+1) 



dr 
2 +(R 



r ui , (2.io) 



with appropriate boundary conditions. These will be 
discussed further in section IIII Bl for the particular case 
of lukewarm black holes. Typically, p^i is the solution 
which is regular at the lower limit of the region under 
consideration (for black holes, the event horizon), while 
q^i is regular at the upper limit of the region (usually 
infinity, but in our case the cosmological horizon). The 
function T w £(r, r') is then given by [15| 



T^(r,r') = C u tPut{r<)qui(r>), 



(2.11) 



where, as usual, r K is the lesser of the two values r, r' and 
r> is the greater. In (|2.11|) . the normalization constant 
C u t is fixed by the normalization condition [l5| 



a 



dqu. 



dr 



dp u 



dr 



1 



2 ,f 



Further properties of the mode functions p^e and q^i will 
be discussed in section ITlIB I 

We now follow standard procedure and choose tempo- 



ral point-splitting, so that 



y and ip = p' . 



Then 7=1 

G e (t,x;t' ,x) 



and, using Pe(l) — 1, we have 
T 



4n 



E - 

n— — oo 

oo 



5^(M + l)CUp^(r) ?a ,/(r), 

(2.13) 



where c 



B. Finite Mode Sums 

Even though the points are separated, the Green's 
function (12.13|) is apparently divergent. This is due to our 
choice of point-splitting and arises in practice because the 
sums over £ in (|2. 13|) do not converge. This problem is 
well-known, arising first in Candelas calculations on 
a Schwarzschild background. This divergence cannot be 
'real' since, by definition, the Green's function must be 
finite when the points are separated. The apparent diver- 
gences are removed by subtracting from (|2.13|) a suitable 
multiple of the delta function (which vanishes when the 
points are separated). The answer is well-known, and 



given by [U[TH: 
G e {t,x;t' ,x) 



T 

47T 



E 



E 

e=o 



(21 + 1) Cuip u i(r)q M t(r) ~r 

rj 2 . 



(2.14) 



Although this answer has been previously calculated, in 
this section we would like to take a different approach 
to deriving the result (]2 . 14[) based on dimensional reduc- 
tion. This new method may prove to be useful in more 
complicated situations (such as Kerr black holes). 

We begin by writing the general Euclidean Green's 
function l|2.7p in the form 

rp oo 

G E (x;x') = j- E ex P N ( T - r ')] 

n— — OO 

xg u (r,6,<p-y, (2.15) 

where 



(2.12) (r, 8, ^ r', 9', <p') = (21 + 1) T^(r, r')P£(cos 7 ). 



(2.16) 

Using the wave equation (|2.5p . the function 
Q u (r, 8, ip; r' ,8', ip') satisfies the differential equation 



S(x, x') 



_L__9 

r 2 dr 



f' 



1 



d 



i d 2 g^ 



dip 2 



r 2 sin 2 i 



r 2 sin 8 88 

2 



m 2 +£R 



89 

Gu>- 
(2.17) 



Eq. (|2.17p looks very much like a wave equation in three 
dimensions, with a potential term: 



V 



f 



m 2 + £R. 



(2.18) 



In fact, by multiplying (|2.17p by / , we obtain the wave 
equation on the three-metric 

d~s 2 = dr 2 + r 2 fd8 2 + r 2 f sin 2 8 dip 2 , (2.19) 

namely: 

V<V*-7(x)] g u (x;x') = -~g-^{x)5 3 (x,x'), (2.20) 

where the covariant derivatives are with respect to the 
three-metric (|2.19[) (acting on x) and g is the determinant 
of the three-metric (|2.19j) . with the potential 



V(x) 



V _ u m 2 +£R 

7~P + 7 ' 



(2.21) 
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where it should be stressed that R is the Ricci scalar of 
the original, four-dimensional, metric (|2.4[) . 

The three-metric (|2.19|) is curious. We emphasize that 
it has no physical significance, and in fact it has a cur- 
vature singularity at an event horizon where / vanishes. 
Furthermore, the potential (I2.21|) is also divergent at a 
horizon. On the other hand, if we are dealing with an 
asymptotically flat metric (|2.2p . then the three-metric 
(|2.19p is also asymptotically flat as / — > 1 at infinity. 
The three-metric should simply be regarded as a useful 
mathematical tool. 

From Eq. (|2.20[) . we can see that Q u (x\ x') is indeed a 
three-dimensional Euclidean Green's function for a scalar 
field on the three- metric (|2. 19|) , but with an unusual po- 
tential (|2.21j) . However, the potential does not affect 
general form of the singularity structure of the Green's 
function, which has the usual Hadamard form in three 
dimensions 129, ugl : 



U(x,x') 
[2a(x,x')]i 



W(x,x'). 



(2.22) 



Here, 2a(x, x 1 ) is the square of the geodesic distance be- 
tween two closely separated points x, x', and U(x,x'), 
W(x,x') are symmetric biscalars which are regular in 
the limit x — ► x', whose precise form will depend on the 
potential V(x). Note that there is no logarithmic term 
in the Hadamard expansion (|2.22j) as we are currently 
working in three rather than four dimensions. The bis- 
calars U(x,x') and W(x,x') can be expanded in terms 
of a(x,x') using standard methods [H, [27j . For our 
purposes here, we only require the lowest order term in 
U(x,x'): 



U(x,x') = l + 0(a), 



(2.23) 



which does not depend on the potential V (|2.2ip . 

We now choose a point-splitting for the three- 
dimensional Green's function Q^. We choose r — r' and 
ip — ip' . Since the metric (|2.19[) is spherically symmetric, 



J 



we can, without loss of generality, fix 0' = 0, and then 
cos 7 = cos 9 (|2.8p and (|2 . 1 6[) takes the form 

OO 

= 1) a#w(r)Mr)P^(cos0). (2.24) 

Although we have brought the radial co-ordinates to- 
gether, this sum is still finite because of the P^(cos#) 
terms. For this point splitting, to leading order we have 



2a = r 2 f9 2 + 0(6> 4 ). 
Therefore the Hadamard form (|2.22p reads 



(2.25) 



Quj (r, 6, ip; r, 0, (p) = — : h finite terms, (2.26) 

so that, for small 0, we have 

y2(2£+l)C u i Pa ,e{r)q u e(r)Pe(co S 6) = — T -+0(l). 

This clearly shows that the sums over £ in (|2.f 3p diverge. 

We find the appropriate subtraction term to render the 
sums over I finite using the identity 0]: 

OO 

Y / Pi(cos6) = - + 0(8). (2.28) 

1=0 

Multiplying this by l/(r/ 1 / 2 ) and subtracting, we there- 
fore obtain: 



E 



{2i + l)C ull p ut {r)q ut {r) 



1 

rfz 



P e (cos9) = 0(1). 



(2.29) 

We may now take the limit 9 — > and find that the sums 
over i in (|2.14p are finite as required. 



C. Renormalized Expectation Value 



From (I2.f 4p . we now have an expression for the unrenormalized expectation value 





\t 


= lim < 






4tt 




It 


= lim < 




1 


2vr 


T 


OO 



n— — oo 

oo 



(2£+l)C u i Pue (r)q ue (r) 



1 

rfi 



i=0 



[2i + 1) C ut Pui(r)q ut (r) r 

rf 2 . 



4?r 



E 

£=0 



(2£ + l)C omi (r)q t{r) 



I 

rfi 



(2.30) 
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To renormalize this expression, we subtract the usual divergent terms (4> 2 }div before taking the limit e — > [2£ 



1 



/div = 



1 



8tt 2 ct 8tt 2 



1 



167T 2 967T 



(2.31) 



where cr Q = <r ;Q , the quantity C is Euler's constant, and R a p is the (four-dimensional) Ricci tensor. If we are 
considering a massive scalar field, then the constant \i is simply equal to to, the mass of the field. However, for a 
massless scalar field, the constant [i is arbitrary [l[ . It corresponds in the massless case to a finite renormalization of 
terms in the gravitational action. 

For our particular point-splitting, the quantities appearing in (|2 . 31[) have already been computed [l], [lj| [28| : 



a 
a" 



(7 V = 0. 



(2.32) 



The subtraction terms (|2.31|) then simplify to [If 



6 2 )div 



1 1 

47r 2 e 2 / + 8^ 
1 d 2 / 



2 V 4 



i 4f 



967r 2 dr 2 487r 2 r dr 



1 

16tt 2 ' 192tt 2 / \dr 



(2.33) 



In order to subtract (|2.33[1 from (|2.30p . we first need to write (|2.33| in terms of mode sums. This is done using the 
following identities, valid for small e and any n > [3,[3]: 



— k 2 ncos(riKe) — — + 0(e 2 ); 
n=l 

cos(n/ce) | Q 2 



(2.34) 



We choose k = 2-7rT, so that u> in Eq. (|2 . 30[) is equal to nn. Then the subtraction terms (|2.33[) become 



"l - 9T / w cos(we) - — j 
Ait 2 f ^ 87T 2 

J n=l 



1 

S^ 2 

o 

m 



4f 



c 



cos(we) 



n=l 



l d 2 / i 4f 



16tt 2 192tt 2 / V*V 96tt 2 dr 2 48?r 2 r dr 48tt 2 / 



+ 0(e 2 ). 



(2.35) 



This is now in a form suitable for subtracting from (|2.30p . and the limit e — > can then be taken to give the final, 
renormalized, expectation value [![: 



)ren — (0 2 ) analytic + (</> 2 ) 



numeric 



(2.36) 
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where 



2\ 

/analytic 



/numeric 



m 2 1 fdf 

16tt 2 ~ 192tt 2 / [dr 



1 

'8^2 



i d 2 f 



OO I oo 



27r E \ E 



{2l+l)C^p^(r)q^(r) 



71=1 (.£=0 L 
oo 



1 


df 


48?r 2 r 


dr 






4k 2 y 


i 


i 







uj 1 
7 + 2^ 



4tt ^ 



(2l+l)C epae(r)qoe(r)- 



V7 



(2.37) 



(2.38) 



As the name suggests, (</> 2 ) analytic has a simple form 
which can be easily computed for any k and metric func- 
tion /. On the other hand, the quantity (</> 2 ) n umcric can, 
in general, only be computed numerically as the mode 
equation (|2.10[) needs to be integrated. Solutions to the 
mode equation are not known in closed form for general 
uj even for Schwarzschild black holes. One might hope 
that ((ft 2 ) analytic is a good approximation to (0 2 ) ren , at 
least in some region of the space-time. This is discussed 
in [l[ for the case of Reissner-Nordstrom black holes, and 
we will examine this issue for lukewarm black holes in 
section IIIIl 



D. Computation of {(ft 2 ) 



numeric 



Before we can compute (</> 2 ) n umeric, further work is 
needed. The mode sums in ((ft 2 ) numeric converge so slowly 
as to render their numerical computation impractical. 
We therefore employ a WKB-like approximation, which 
will give the large uj, I behaviour of the mode sums in 
(I2.38p . We have found it simplest to use the WKB ap- 
proach of Howard [3j , which is different from that used 
in Ref. [JOjl. 

Define a new function fi u i(r) for u> ^ by 



/3 u e(r) = CuePwiirfqutir). 



(2.39) 



Then p u g satisfies the following differential equation (cf. 
0): 



Pull 



1 



%Xu>t 



f 



1 d 2 (V737i) 



V 



xle \ VP7e d( 2 

(2.40) 

where we have defined a new independent variable C by 



d _ 2 d 



(2.41) 



and the functions Xuil{ r ) and v( r ) are given by 
Xulix) = \ uj 2 r 4 + (i + 



r 2 f 



1 



rj(r) = --fr 2 + (m 2 + £R) fr 4 . (2.42) 

We are looking for the behaviour of fiutir) when either 
uj or I are large, that is, when x^i( r ) is large. This is 
found by inserting a fictitious parameter e into (|2.40p as 
follows: 

1 d 2 (v%Tt) 



1 

2x~ 



1 - 



1 



V 



(2.43) 

and then expanding U £ in inverse powers of s: 

Mr) = Mr)+e- 2 f3 luje (r)+e~ i p2u>i{r) + .. . , (2.44) 

finally setting e = 1 at the end of the calculation. 

It is straightforward to read off from (|2.43[) that Pa^tir) 
has the simple form 

Mr) = ^TT- (2-45) 

The next term in the expansion Piu,i(r) is more compli- 
cated but not difficult to calculate: 

Mr) = -^Jjr ( A ixtt + B ^li + Ci) , (2.46) 



where 
Ai 



F2 r 4 



64 X ; 

-fr 



2 I 

dr 



12 fr^f + 4/ 3 - 4/ 2 



-16r 2 £R/ 2 + Af 



dr 2 



dr 

+ 16r 2 m 2 / 2 



-16/ 2 r? + 16/ 3 



dr 



4,*+rMf 



4/r 



(2.47) 
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It can be seen from (12.4512. 46[) that, for large XojI ( or , to be considered separately. When u> = 0, from (|2.42[) we 



equivalently, for large uj or i if / ^ 0), 



have xoi = + 5) an d therefore (|2.45I2.46[) become 



(2.48) 



and we find similarly that fii^t ~ xJi e ^ c - If we consider 
the sum 



00 00 



Axw(r) 



/8io/(r) 



2r/* 



5] ^ (2* + 1) [CUp^(r)^(r) - /3 ^(r) - /3 



64/J 



(2.50) 



n=l £=0 



(2.49) Therefore (2£+l)(3 oi{ r ) = (rf*) * and this P art ° f thc 



WKB approximation has already been subtracted from 



then the summand is 0(£ 4 ) for large I. This is the 

behaviour we observe when we compute the summand theTum" Yn thi7case we consider 
numerically in section IIII Bl When this is summed over 
t, we obtain a summand which is 0(ui~ 3 ) and therefore 
converges rapidly. Therefore it is sufficient to subtract 
just Pout an d fliul from the mode sums in 1 



(2* + 1) [C uPu (r)q M {r) - (3 Qoe (r) - (3 10i (r)] 



/numeric 



£=0 



(2.51) 



The contribution to (|2.38|) from the u — modes needs and again the summand is 0(£ 4 ) for large I as required. 
I 



We therefore write 



/ numeric 



T 
27 



OO f OO 



E E(^ +1 ) v»& 

ut.(r)quje(r) - (3 0ui i(r) - ft u <(f 



E 

.Hi 

4tt 



(2* + l)[^(r)+/Wr)]-4r 



w 1 

7 + 2ZJ 



E ( 2 * + J ) [CoiPM(r)««(r) - Axw W - 1M (r)] + E ( 2 ^ + 1) /M»0 I , (2.52) 



«=0 



where we have used the fact that (2£ + l)Pooe(r) — (rf*^ . The first and last lines in (|2.52[1 are now amenable to 
numerical computation. The final sum is readily computed using (|2.50|) to be 



°° 3 00 / 1 \ ~ 2 23 

E (2* + 1) Mr) = -^-tM E 2 K + o = -tHt^ 1 



€=0 



64/ 



£=0 



64/3 



(2.53) 



Since we have analytic expressions for Pouji( r ) an d Piuje(r), we next examine the term 



UJ 1 

/ + 2oZ 



m + [ ( 



6 



i? 



(2.54) 



The sums over € are most easily found using the Watson-Sommerfeld formula (see, for example, [4]), valid for any 
function analytic in the right-hand half plane: 



f 1 1 \ f 2 / I 



(2.55) 



Using (|2.55[) . we write, for w 7^ 0, 



E(2^ + l) 



A)w£(r) TT 

rf*. 



Io(w,r) + Jo(w,r); 



J2(2£+l)(3 luje{r) = 7 x (a;,r) + Ji(w,r); 



(2.56) 
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where the 7j(o;,r) are the integrals from the first term in (|2.55[) and the Ji(ui,r) are the integrals from the second 
term in (|2.55p . The I^(uj,r) integrals are easily calculated for w ^ (here I = A — 1/2): 



Substituting these results in (|2.52[) gives 



2X(3 ^(r) - —r 
2A/3 w (r) d\ - 



1 

2uj 



R 



24r 2 o 



(2.57) 



/ numeric 



T 
2^ 



oo f oo 



(uj, r) + J\{w, r) 



1 



n=l kfcO 
T 



24r 2 w 



. 1=0 



64/ 



(2.58) 



Next we examine the Ji(w,r) integrals (here £ — iX — 1/2): 



J (w,r) = -K{i 



1 + e 



2ttA 



2Az/W(r) r 

rf 2 . 



dAL Ji(w,r) = » 



4A 



1 + e 



2ttA ' 



Pi^{r)d\\. (2.59) 



We consider each of these in turn. Firstly, for Jo(w,r), substituting in the form of 0o u (.{r) (|2.45[) gives 



J ( w ,r) = — T & 
where we have defined 



2? 



o 1 + e 



2ttA 



dA 



2A 



o (1 + e 27rX ) (a 2 — A 2 ) 



■dA 



2A 



a (l + e 27rA )(a 2 -A 2 )- 



-dA 



f* 



(2.60) 



(2.61) 



The first and third integrals in (|2.60p do not contribute, leaving just the second integral. The integrand is integrable 
but not regular at A = a. Furthermore, this integral is not known in closed form and will need to be computed 
numerically. It is easier numerically to have a regular integrand, so we integrate once by parts to obtain 



J {u),r) = — - 



For large w, this integral behaves like 



/ f Jo 

Jo(w,r) 



_ A*\ 2 
1_ a 2 



2ttA 



(l + e 27rA ) 



■dA. 



1 



24wr 2 



+ 0{u;- 3 ). 



We therefore take the — (24wr 2 ) 1 term from (|2.58|) and consider instead 



Jo(w,r) 



1 



47T0J 



/ 24r 2 w / J Q 



A 2 V 



= 2ttA 



( 1+e 27rA) 



dA - 0(uj- 3 ) 



as ui — > oo. This means that the sum X^^Li Jo(w,r) can be computed separately from the other terms in 
The integral Ji(w,r) is more complicated. Using the form (|2.46[) of /3i w ^(r), we find: 



Ji(u,r) = - 



1 



16r/i 



r 4 Aii^i(w,r) + \B 2 K 2 (u,r) + -^dK 3 (uj,r), 



(2.62) 



(2.63) 



(2.64) 



/numeric 



(2.65) 



where we have defined new integrals 



-5R 



K 2 (tu,r) = -SR 
a 



(1 - g 2 ) 2 (1 + e 27rog) 

9 

o (1 - q 2 )^ (1 + e 27 ™?) 

q 

o (1 - q 2 y- (1 + e 2 ™i) 



dq 
dq 
dq 



(2.66) 
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and q — A/a. The integrands in all three integrals have branch points at q = ±1, so we cut the plane along the 
interval [—1, 1] and consider the contour shown in Fig. [TJ 




FIG. 1: Contour used for computing the integrals Ki(uj,r) (|2.66|) . 



The contribution to each Ki(u>, r) from the contour 72 is zero, so we write, for each i = 1, 2, 3: 

Ki(u,r) = - lim [L^oj, r) + M^w, r)] 

a e^O 



(2.67) 



where each Li(uj,r) is the contribution from the contour 71 and Mi(u>, r) is the contribution from the contour j e in 
Fig. [TJ We illustrate the procedure for calculating these by considering L\(uj,r) and M\(u>, r). The method works 
similarly for K 2(1-0, r) and K^(lo, r) but is more complicated. 
Considering Li(u,r) first, we have 



Li(w,r) = H 
Integrating by parts gives 

Li(uJ,r) 



■ (1 - g 2 ) 2 (! + e 27ra 9 ) 



■ dg 



(1-g 2 ) 2 (l + e 27rQ 9) 



■ dg. 



(2.68) 



1 1 „ 

— 1=-, z — ; h 27ra 

2Vi(l + e 2 - Q ) 2 Jo 



l-E 



1 - g 2 ) 2 e 2 ^ 
(l + e 2 ™ 1 ?) 



• dg + 0(ei\ 



(2.69) 



Note that the first term in (|2.69|) is divergent as e — ► 0, and that we have now isolated a finite integral (although the 
integrand is not regular). We next turn to M\(oj, r): 



Mi(w,r) = 9? 



1. (1-g 2 ) 2 (l + e 2 ™?) 



■ dq 



Along 7 e , we have 

q - 1 = -ee i,? , < < tt. 
Changing the variable of integration in M\(ui,r) to J?, and computing the integral, we find 

1 



2Ve(l + e 27rQ ) 



■0(e* 



(2.70) 



(2.71) 



(2.72) 



Therefore, adding Li(o;,r) and Mi(uj,r) will give a finite quantity and we can take the limit e — > 0. It is helpful for 
numerical computation to integrate by parts again to give a regular integrand: 



1 



7r 2 e 2 ™ 



2a ' (l + e 27ra) 2 

where we have returned to the original variable A = aq 



47T Sill 1 1 > 



1 



(l + e 2 ^) 



dA, 



(2.73) 
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The calculation proceeds similarly for K^ui, r) and K^{uj,r), except that more integrations by parts are required. 
In each case, the divergences in Li(u/,r) and Mj(w, r) cancel to give the following finite quantities: 

1_ 8^a r a / _)?\* e 2 " A (-l-e 4 " A + 4e 2 " A ) 
6a + 3 Jo ^ aV (1 + e 2 - A ) 4 



^,r) = -- + —/ ^— — 3 ' > dX; (2.74) 



1 167r 4 a [ a ( X 2 Y e 27rA (-l + lle 27rA -lle 47rA + e 67rA 



K 3 (a;,r) = / A 1 5- ^ e L d\. (2.75) 

V ' 10a 15 J \ a 2 J (1 + e 2 - A ) 5 V ' 

It is straightforward to verify that, for each i, the integrals Ki(w,r) are 0(w~ 3 ) for large u, so that the sums over n 
required in (|2.58j) will converge rapidly. 

At this stage it is helpful to bring together our results for (</> 2 ) numeric- From (|2.58|) we have 



^ ) numeric 



S + — Y,[M",r) + Ji(u>,r)], (2.76) 



2tt 

n=l 



where 



rp OO f CO | 

S = 2~ E 1 E ( 2 * + !) [^£P^(r)o^(r) - /3 ^(r) - /3 lt ^(r)] 

+ i E ( 2£ + !) [GxPM(r)gw(r) - /3 m (r) - /W)] - ^-3-^ , (2.77) 

and the numerical integrals J^{uj,r) and J\{uj,r) are given by (|2.64p and (|2.65)1 respectively (with the Ki(tv,r) given 
by (|2.73H2.75"|) ). In the following we shall refer to S as the 'mode sum'. 



III. NUMERICAL RESULTS FOR LUKEWARM BLACK HOLES 

We now turn to the computation of (</> 2 ) r0 n for a specific example, namely lukewarm black holes [2(|. We briefly 
review the key features of these black holes before describing the results of our calculation. For many of our numerical 
computations we have used standard routines [2^ | . Therefore we omit much of the detail of the numerical methods 
used, which are comprehensively discussed in Ref. [30] - 



A. Lukewarm Black Holes 

Lukewarm black holes are a particular type of Reissner- 
Nordstrom-de Sitter space-time, with metric given by 
(|2.2|2.3p with M = Q. For AM < ^/Jfk, there are three 
distinct horizons (a black hole event horizon at r = r +1 
an inner (Cauchy) horizon at r — r_ and a cosmological 
horizon at r = r c ), given by 



r_ = 


L 
2 


-n 


1 AM 

-v 1+ — 


r+ = 


L 
2 


1 - 


/ AM 




r c = 


L 
2 


1 + 


/ AM 





where 

L = \Hr, AM < L. (3.2) 



The Penrose diagram for this space-time can be found in 




In this case the event and cosmological horizons have 
equal surface gravities: 



«+ = Ke = -^y 1 - — , (3.3) 

which means that they have the same temperature T = 
n + /2tt. Here we are interested in the region between 
the event and cosmological horizons, which has a regu- 
lar Euclidean section, with topology S 2 x S 2 19]. Our 
computations are for a thermal state at the natural tem- 
perature T. 

In the following sections, we show plots of the con- 
stituent parts of (</> 2 )ren for the specific case of M = Q = 
0.1-L, although the main features are the same for other 
values of M. For the remainder of this section, all di- 
mensionful quantities (r, M, etc) are given in units of 
L. 
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Mode Sum 



We begin our numerical analysis by calculating the 
mode sum (I2.77[) . Firstly we need to find the modes 
themselves, by integrating the mode equation (|2.10[) 
which is satisfied by p u i{r) and (^(r). This is done using 
standard shooting techniques [29| • 

The mode equation (I2.10|) has regular singular points 
at the event and cosmological horizons. Using the stan- 
dard Frobenius method, we write the mode functions as 
power series: 



(3.4) 




o.o 
Event 
Horizon 



0.2 



0.4 

(r-r + )/L 



Cosmological 
Horizon 



where S^e is either p u i(r) or q u e(r), and x = r — r + > if 
we are considering the behaviour near the event horizon 
(x = r c — r > if we are considering the behaviour near 
the cosmological horizon) . By looking at the lowest order 
term in (|2.10[) we obtain the indicial equation [l5l ]: 



(3.5) 



where ro is either r+ or r c as applicable. In our mode 
sum (|2.77p . we are only summing over modes for which 



lu = 2imT = tik+ = 



dx 



ra = 0,l,2,. 



(3.6) 

Therefore v = ±n/2, and the roots of the indicial equa- 
tion differ by an integer (n > 0) or zero (n = 0). There- 
fore the linearly independent solutions of the mode equa- 
tion ([2T0]) are: 



Si u t = ^2cnx l+ 2] 



i=0 



S2ut = * + 1C u e.Si W i\n 



i=0 



(3.7) 



where IC^i is a constant which definitely does not vanish 
if n = but may possibly vanish if n > 0, and ro = r + 
or r c as applicable. The mode functions are chosen so 
that Put is regular (that is, has the form Si^e) near the 
event horizon, while q^g is regular near the cosmological 
horizon. From (|3.7p , it is clear that q u t will diverge (that 
is, have the form S^ui) near the event horizon, and p u e 
will diverge near the cosmological horizon. This can be 
seen in the examples of mode functions plotted in Fig. [5] 

For n > 0, it is clear that the mode functions p w f(r) 
will vanish at the event horizon and that q u e(r) will van- 
ish at the cosmological horizon. However, for n — 0, we 
have that pot(r) is regular but non-zero at the event hori- 
zon (and similarly for qoi{r) at the cosmological horizon). 
In this case poe(r) will still diverge at the cosmological 



FIG. 2: Examples of mode solutions of the equation (|2.10|l . 
with ui = 2-7rT and I = 0. The function p u e(r) vanishes at 
the event horizon (at the left-hand-edge of the plot) and di- 
verges at the cosmological horizon (at the right-hand-edge of 
the plot), while the mode function q^i(r) vanishes at the cos- 
mological horizon and diverges at the event horizon. For these 
mode functions, we have taken the scalar field to be massless 
and conformally coupled, and the black hole metric (|2.2l2.3p 
has parameter M = Q = 0.1L. The radial co-ordinate on the 
horizontal axis is in units of L. 



horizon because of the ln(x/ro) term in S2ut l|3.7|) (note 
that TCqi cannot be zero because there must be two lin- 
early independent solutions of the mode equation). It is 
straightforward to show, from the mode equation (|2.10p . 
that the mode functions are monotonic and do not have 
zeros except possibly at a horizon. 

To find the jw(r) mode functions, we start the numer- 
ical integration of the mode equation (|2.10[) just outside 
the event horizon, using as many terms in the power se- 
ries expansion ()3.7H of Si^e as required for the desired 
accuracy. As in [1 51 ] . there is a complicated recurrence 
relation (involving seven terms in general) for the coef- 
ficients in the power series expansion (|3.7j) . which we do 
not reproduce here. Starting with ao = 1, this recur- 
rence relation is used to compute the power series expan- 
sion. We then integrate outwards towards the cosmo- 
logical horizon. For the qu>t{r) mode functions, we start 
integrating just inside the cosmological horizon, and in- 
tegrate towards the event horizon. We use the Bulirsch- 
Stoer method [29[ of integrating the differential equation 
P-lOP because the mode functions are not oscillating and 
because of the high degree of accuracy required. 

Looking at the mode sum (|2.77p . it can be seen that 
for large u> or the rapid convergence of the sum is 
dependent on subtracting very nearly equal quantities. 
Therefore it is imperative to calculate the mode func- 
tions Put(r), quii(f) with great accuracy. We there- 
fore used quadruple precision throughout our calcula- 
tions. The normalization constant C u i is computed from 
(I2.12p . The constancy of the C^i as calculated from 
(I2.12p for different values of r represents a good check 
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on the accuracy of our results. For the particular case 
of Q = M = 0.1L, we find that C^g remains constant 
to within 1CT 26 \C ul \ for < n < 25 and < £ < 600, 
which gives us enough modes to get good convergence of 
the mode sum (|2.77[) . 

Once we have the mode functions Puig{r), q u £(r), as 
the WKB approximants Aw(r) (|2~45|) and /?w(r) (f2T46|) 
are known analytically, we are able to find the mode sum 
(|2.77p . We found that the summand was 0(£~ 4 ) for large 
£ as predicted in section III Dt We found that the sums 
converged more quickly if we summed over n first and 
then £. Convergence of the sums was speeded up by using 
a Shanks transformation (see appendix D of j3j). The dis- 
advantage of the Shanks transformation is that one loses 
accuracy (typically the accuracy after the transformation 
is half that before) , but we are nonetheless confident that 
our final answers are accurate to five significant figures. 
The form of the mode sum (|2 . 7T[) will be shown explicitly 
as a function of r in section MI Dl 



C. Numerical Integrals 

The remainder of (</> 2 ) numeric involves sums of integrals 
which have to be computed numerically (|2.76p : 



0.3 




o.o 0.2 0.4 0.6 Cosmological 

Event (r-rJ/L Horizon 

Horizon 



FIG. 3: Sums over n of the Ki(uj, r) integrals (|2.73H2.75|) as 
functions of r. Note that these are independent of the mass 
or coupling of the scalar field. The radial co-ordinate on the 
horizontal axis is in units of L, and we consider the black 
hole metric (|2, 212.3ft with parameter M = Q = 0.1L. It can 
be seen that all three sums vanish at the event horizon (left- 
hand-edge of the plot) and the cosmological horizon (right- 
hand-edge of the plot). This is in agreement with analytic 
work to be presented in section llV Al 



D. Final Answers 



(tumeric = W~ £ IM^) + Ji(w,r)] , (3.8) 

where Jo(w,r) is given by (|2.64p and Ji(uj,r) is given 
in (|2.65p as a combination of three integrals, Ki(ui,r) 
(I2T731 . K 2 (u,r) (571 and£ 3 (w,r) (12~T5D . Each of these 
integrals depends on a (|2.6ip . which in turn depends on 
u> = 27mT, for n — 1, . . . (note that the numerical in- 
tegrals are only relevant for u ^ 0, as the modes for 
u! = were considered separately in section Hi D|) . For 
each value of n, the integrals are straightforward to com- 
pute, and then summed over as many n as required for 
convergence. The summands are all 0(n~ 3 ) as n — > oo 
so convergence is relatively rapid. As an example of our 
results, we plot in Fig. [3] 

oo 

£tfi(w,r) (3.9) 

71=1 

for i = l,2, 3, showing how these sums vary with r. 

From (12.73H2.75f . these sums are independent of the 
mass or coupling of the scalar field, but do depend on the 
metric of the black hole through the function / (|2.3p . It 
can be clearly seen in Fig. [3] that all three sums vanish 
at both the event and cosmological horizons. We will 
show analytically in section IIV Al that this is always the 
case, independent of the temperature of the thermal state 
under consideration. 



The final answers for (0 2 ) re n are now computed by 
combining the mode sum (|2.77p with the numerical in- 
tegral contribution (|3.8p to give (</> 2 ) numeric (I2.38p . and 
then adding this to (0 2 ) a nai y tic (|2.37[) . 

We begin this section by examining the contributions 
to (0 2 )numeric, namely the mode sum (|2.77[) . and the two 
sums over numerical integrals 

T t 

-5>(.,r), ^£<M^h (3-10) 

n—l n—1 

These are shown in Figs. HH7] for a lukewarm black hole 
with M = Q = 0.1L. The quantity labeled 'Total' in 

FigS. SHE is ) numeric- 

In Figs. [4] and [5] we consider a massless, conformally 
coupled, scalar field, while in Figs. [5] and [7] the field is 
still conformally coupled but has mass m — L. We have 
chosen a large value of the mass so that the differences 
between this and the massless case are very clear. 

Examining first the massless case, from Fig. 0]we can 
see that all three contributions to (</> 2 )numcric have peaks 
near the event horizon of the black hole, and that all 
three parts become very small at the cosmological hori- 
zon. The behaviour near the event horizon is clearer in 
the close up in Fig. [5] From Fig. [S]it can be seen that 
the contributions coming from sums over the numerical 
integrals p.lOp tend to zero as the horizon is approached, 
and the dominant contribution to (<^ 2 ) n umeric comes from 
the mode sum S (|2.77p . Furthermore, this appears to re- 
main finite as the horizon is approached. 
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0.03 
0.02 
0.01 
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M & ter) 
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^(mfi 


Total 



0.0 

Event 
Horizon 



0.2 



0.4 

(r-r + )/L 



06 Cosmological 
Horizon 



0.02 



0.00 



-0.02 



-0.04 




0.0 

Event 
Horizon 



0.4 

(r-r + )/L 



06 Cosmological 
Horizon 



FIG. 4: Contributions to (</> 2 } n umcric (|2.38f) from the mode 
sum E (|2.77|l and the sums over the numerical integrals (|3.10|l , 
as functions of r. Here, the quantum scalar field is massless 
and conformally coupled. The radial co-ordinate on the hori- 
zontal axis is in units of L, and the black hole metric (|2.2I2.3|) 
has parameter M = Q = 0.1L. All quantities appear to be 
regular near the cosmological horizon. Their behaviour near 
the event horizon is shown in close-up in Fig. [S] 



0.03 
0.02 
0.01 
0.00 
-0.01 
-0.02 
-0.03 




E 

£ n Ji(». r ) 
Total 



-0 04 



0.000 0.001 0.002 0.003 0.004 0.005 0.006 

Event (r-rJ/L 
Horizon 



FIG. 5: Close-up, near the event horizon, of the contributions 
(shown in Fig. [3]) to (</f> 2 } n umcric (|2.38p from the mode sum 
E (|2.77[l and the sums over the numerical integrals (|3. lOfl . 
when the quantum scalar field is massless and conformally 
coupled. It is clear that the contributions from the sums over 
the numerical integrals (|3.10[l vanish at the event horizon. 
Near the horizon the mode sum E (|2.77p is the dominant 
contribution and this appears to remain finite as the horizon 
is approached. 



When the scalar field is no longer massless, we see the 
difference in behaviour clearly in Fig. [5] Near the cos- 
mological horizon, the contributions to (</> 2 ) numeric from 
the sums over the numerical integrals (|3.10[) still vanish 
as the horizon is approached, but now the mode sum E 
(|2.77p (and therefore (0 2 ) n umcric) diverges. The close-up 
near the event horizon in Fig. [7] confirms that this also 
happens near the black hole event horizon. 

We will show in section IIVI that the contributions to 



FIG. 6: Contributions to (<^ 2 ) nU moric (|2.38[) from the mode 
sum E (|2.77|l and the sums over the numerical integrals (|3.10|l , 
as functions of r. Here, the quantum scalar field is conformally 
coupled but has mass m = L. The radial co-ordinate on the 
horizontal axis is in units of L, and the black hole metric 
(|2.2I2.3|I has parameter M = Q = 0.1 L. The contributions 
from the sums over the numerical integrals appear to be regu- 
lar near the cosmological horizon, but the mode sum E (I2.77P 
is diverging as the cosmological horizon is approached. The 
behaviour of these quantities near the event horizon is shown 
in close-up in Fig. [7] 



o.oo 



-0.04 - 



Total 



0.000 0.001 

Event 
Horizon 



0.002 0.003 0.004 0.005 

(r-rJ/L 



FIG. 7: Close-up, near the event horizon, of the contributions 
(shown in Fig. |S} to ((f> 2 ) numeric (|2.38p from the mode sum E 
(|2.77p and the sums over the numerical integrals (|3. lOf) , when 
the quantum scalar field is conformally coupled and has mass 
m = L. It is clear that the contributions from the sums over 
the numerical integrals (|3.10[l vanish at the event horizon. 
Near the horizon the mode sum E (|2.77p is the dominant 
contribution and this appears to be diverging as the horizon 
is approached. 



(<fi 2 ) numeric from the sums over numerical integrals always 
vanish at the horizons, whatever the mass or coupling of 
the scalar field, so that the mode sum £ is indeed the 
dominant contribution to (</> 2 ) n umeric near the horizons. 
Furthermore, we will show that S diverges near a hori- 
zon unless the quantum scalar field is both massless and 
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conformally coupled. 

We now turn to the final result, namely (4> 2 ) rC n- In 

FigS. [8H9]we plot {4> 2 ) ana l y tic (0 2 ) numeric 

and their sum (</> 2 ) ren for a massless, conformally coupled 
scalar field, and a conformally coupled scalar field with 
mass to = L. 



~ 0.04 




0.4 

(r-r + )/L 



0.6 Cosmological 
Horizon 



The vacuum polarization 

a2\ 



and its components 

2\ 



FIG 

\*P /analytic 

([OTJl and (0 2 ) 

numeric 

([2738]) (all in units of L' 2 ). 
Here, the quantum scalar field is massless and conformally 
coupled. The radial co-ordinate on the horizontal axis is in 
units of L, and the black hole metric (|2.2l2.3f) has parameter 
M = Q — 0.1L. All quantities appear to be finite near both 
the event and cosmological horizons. Near the event and cos- 
mological horizons, (</> 2 ) analytic is the dominant contribution 
to (c/> 2 ) rcn . 
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FIG. 9: The vacuum polarization {(f} 2 ) reu and its components 
\<P /analytic (|2~37|) and (cf) 2 ) numeric ([2738]) (all in units of L~ 2 ). 



Here, the quantum scalar field is conformally coupled and has 
mass m = L. The radial co-ordinate on the horizontal axis 
is in units of L, and the black hole metric (|2.2l2.3p has pa- 
rameter M — Q = 0.1L. It can be clearly seen that both 
{(p 2 ) analytic and (<j> 2 ) numeric diverge near the event and cosmo- 
logical horizons, but with opposite signs. 



When the scalar field is massless and conformally cou- 
pled, from Fig. [5] it is clear that (0 2 ) a nai y tic, (0 2 ) numeric 



analytic has the same 
(<^ 2 )ren but does not 



and (</> 2 ) r cn are regular at both the event and cosmo- 
logical horizons. The total (4> 2 ) YBn is greatest near the 
event horizon, and slightly negative near the cosmolog- 
ical horizon. The contribution (</> 2 
qualitative behaviour as the total 
represent more than a qualitative approximation to the 
exact quantity. 

When the quantum scalar field is not massless, from 
Fig. [9] we can see that both (0 2 ) analytic and (0 2 ) numeric 
diverge as either the event or cosmological horizons are 
approached, albeit with opposite signs. From our numer- 
ical results, it looks like these divergences cancel to give a 
finite (0 2 )rcn- However, finiteness of this quantity cannot 
be proven by numerical calculations alone. Therefore, in 
the next section, we study (0 2 ) re n analytically. 



IV. REGULARITY ON THE HORIZONS 

Our numerical results, computed in the previous sec- 
tion, indicate that (t4 2 ) rc „ is finite on both the event and 
cosmological horizons, for the thermal state at the same 
temperature as both horizons. In this section we shall 
show analytically that this is indeed the case. We shall 
consider a general, non-extremal horizon at r = rrj, and 
show that (0 2 )ren is regular at that horizon for a thermal 
state at the same temperature as the horizon. 

A. Analytic Part and Numerical Integrals 

We consider first the behaviour of (</> 2 ) analytic at a hori- 
zon r = rrj . In order to analyze the regularity at a horizon 
of (</> 2 ) analytic P-37|) . we split it into two parts: 

/ j2\ 



/analytic 



where 



^ ) analytic 



h 2 \ H 

l analytic 



h 2 ) 1 

/ analytic 



l 



I analytic ) 



(4.1) 



16tt 2 



I 



1927T 2 / 



1 d 2 f 
96tt 2 dr 2 



48n 2 r dr 



1 

'8^2 
C 



48tt 2 /' 
1 



R 



In 



An 2 ) 



(4.2) 



At first sight it looks like (</ ,2 ) an aiytic ^ s divergent at a 
horizon where / = 0, but in fact it is straightforward to 
show that if k is chosen to be ^^f- at the horizon, then 

2 dr 

(fi 2 ) analytic ^ s regular there. This is in accordance with 
expectations: if the temperature T — 2irn of the state 
matches the natural temperature of the horizon we are 
considering, then (</> 2 ) ana i y tic i s regular at that horizon. 
In particular, for the lukewarm black holes, we have that 
(fi 2 ) analytic ^ s regular at both the event and cosmological 
horizons. 
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On the other hand, it is clear that (0 2 ) analytic nas a 
logarithmic divergence as / — ► unless both m 2 = and 
£ = 1/6, that is, the field is massless and conformally 
coupled. We will show in the next subsection that this 
divergence cancels with a corresponding divergence in the 
mode sum E (|2.77|) . to give an overall (<^ 2 ) re n which is 
regular. 

Before that, we examine next the other contribution 



to 



/ numeric ) 



namely 



rp °o 

/^numeric = 7^ [M^ , r) + J^UJ , 

n=l 



r)}, (4.3) 



where Jo(uj,r) and Ji(oj,r) are given, respectively, by 
(|2.64l2.65p . with Ji(uj, r) in turn determined by the inte- 
grals in (|2. 731 1^7751) . It is straightforward to show that, 
for fixed lo and small /, the integrals Jo(uj,r), Ki(u,r), 
K 2 (u,r) and K 3 (cu,r) behave as follows: 



J (w,r) 
K x {u,r) 



7/ 



1920w 3 r4 



48w 3 r 3 
48w 3 r 3 

f 2 



0(f); 

21/j 
3840^ 

7/8 
768w 5 rg H 

49/* 
3840w 5 r^ 



0(fi 



+ 0(/»); (4.4) 



where r — ro is the location of the horizon where / = 0. 
We should comment that the first terms in the expansions 
of Ki{u, r), K2(uj, r) and K^{oj, r) do have r and not r$ 
in the denominator. Therefore they include some sec- 
ond order terms as well as leading order terms. Writing 
the expansions in this way simply makes the expressions 
shorter. From (|4.4[) . it is clear that the sums over the 
Ki((J,r) should vanish at the horizon, regardless of the 
temperature T of the state under consideration or the 
mass or coupling of the scalar field. This is in agreement 
with our numerical results in Fig. 

Combining the Ki(tv,r) integrals to form Ji(cu,r) us- 
ing (|2.65p . we find that the leading order contributions 
from the Ki(u,r) integrals cancel and Ji(w, r) has the 
following behaviour: 



Ji(w, 



7/ 



7680w 5 7 



o(/ 2 ). 



(4.5) 



Recalling that w = 2nirT, n ^ 0, we therefore see that 
both Jo(w, r) and J\{oj, r) tend to zero uniformly in n as 
/ — ► 0, and therefore the contributions to (0 2 ) n umeric 
from the numerical integrals (|4.3[) vanish at both the 
event and cosmological horizons. This is independent of 
the (non-zero) temperature of the state under consider- 
ation, and the mass and curvature coupling of the scalar 
field, and is in agreement with our numerical results in 
Figs. EH3 



Therefore the only parts of (<fi ) ron which need further 
investigation are ((f> 2 ) ^ alytic (|4.2[) . which has a logarith- 
mic divergence as r — ► ro, and the mode sum (|2.77|) 
contribution to (</> 2 ) numeric- If (0 2 )ren is to be regular 
at the horizon, it must be the case that the mode sum 
E (|2.77p also has a logarithmic divergence, and, further- 
more, this must exactly cancel the logarithmic divergence 
in (0 2 ) analytic- We next establish that this is the case us- 
ing a careful analysis of the mode functions. 



B. Mode Sums 



The behaviour of the mode sum E (|2.77[) near a hori- 
zon is by far the most difficult part of (</> 2 ) r0 n to analyze. 
The reason for this is that the sums over t are not uni- 
formly convergent: more and more £'s are needed before 
convergence is reached as we get closer and closer to a 
horizon. This is due to the fact that the WKB approx- 
imants floul(r) (|2.45[) and /?w(r) (I2.46|) do not depend 
on I when / = 0, as in this case Xuii(r) Ij2.42p reduces to 



Xwt{r) 



(4.6) 



Furthermore, the WKB approximants do not correctly 
reproduce the ln[(r — ro)/Vo] behaviour seen in the field 
modes Q3.7P near the horizon. 

Suppose we are considering a horizon on which the 
mode functions p u i{r) are regular and the q u e(r) diverge. 
The analysis below will be equally valid if we are on a 
horizon where the 2w( r ) diverge and the q^ii?) are reg- 
ular, but with the p's and g's swapped over. Then near 
the horizon, the mode functions p u t{r) and q u e(r) will 
have the following behaviour (|3.7p : 



(r - ro)* + 0(r - r ) 1+2 ; 
(r - r )~% + 0{r - r ) 1_ * 

+AC w <(r-r )*ln' ' " ''" 



+o 



(r - r ) 



1+i 



hi 



ro 
r - r 
ro 



(4.7) 



for some constant K. u t, where w 
the horizon, we have 



2nnT. Therefore, near 



Cuip^i^q^r) ~ CuilCuiir - r ) n In 
+finitc terms. 



r - r 
ro 



(4.8) 



From (|4.8p . it is clear that the only contribution to the 
mode sum E which is important near the horizon is that 
for n = = to, since this contains a ln[(r — ro)/ro\ 
term which diverges at the horizon. For n > 0, the 
(r — ro)" ln[(r — ro)/Vo] term is subleading compared to 
the finite terms. The finite terms in (|4.8p for n > are 
described well by the conventional WKB approximation 
described in section |TTD] and their contribution to the to- 
tal mode sum is finite as the horizon is approached [23|. 
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To analyze the behaviour of the mode sum E, what 
is therefore required is an approximation for the mode 
functions with n — uj — 0, near the horizon, which is 
uniformly valid in I. Such an approximation was first 
developed by Candelas [|[ for Schwarzschild black holes, 
and subsequently extended in Ref. [23[ for the Reissner- 
Nordstrom black hole. Here we generalize the approach 
of Ref. [23| for general metric function f(r). We focus 
on the case of a horizon at r = To, where poe(r) is regular 
and qoe(r) diverges, but the analysis proceeds similarly 
for a horizon where q$i (r) is regular and p$i (r) diverges. 

We begin by writing the mode functions in terms of 
modified Bessel functions [2311 : 



Poe(r) 

m{r) 



r(r) 
r(r) 



(4.9) 



the mode functions near the horizon. These are then, in 
turn, substituted into the mode equation (|2.10[) . which is 
then expanded in powers of x. The lowest order term in 
(|2.10[) gives the following expression for the constant B: 



B 



(4.14) 



where Rq = R(ro) is the value of the Ricci scalar on the 
horizon. The expression (|4.14|) generalizes that found in 
[23| for the Reissner-Nordstrom case. As in (23|, the next- 
to-leading order term in the mode equation (|2.10p gives 
a complicated expression involving y\ and y 2 - Using the 
fact that both y\ and y 2 must remain bounded as t — > 00, 
we obtain 



where T(r) and O(r) are functions which we shall define 
very shortly. The mode functions (|4.9[) are required to 
satisfy the differential equation (|2. 10|) and the normaliza- 
tion condition (|2.12|) . Considering the normalization con- 
dition (|2.12p first, this is identically satisfied with Coe = 1 
if r(r) is defined to be 



T(r) 



»(r') 
fir') 



dr' 



(4.10) 



whatever the (as yet undetermined) function Q(r). With 
this definition of T(r), substituting (|4.9[) into (|2. 10[) gives 
a differential equation for £l(r) which can be found in 
23], the function fl(r) is defined to be 



[23J. Following 



(4.11) 



where y(r) is a function which is regular at the horizon. 
To analyze the behaviour of the mode functions near the 
horizon, we expand f(r) and y(r) near the horizon as 
follows: 



f(r) - x[f 1+ xf 2 + 0(x 2 )] ; 

y(r) = B[l + x yi +x 2 y 2 + 0(x 3 )] ; 



(4.12) 



yi 



0{t 



00, 



(4.15) 



which means that, as i — > 00, we have 



1 



B l = i {I + 1) + - + m 2 rl + ( £ - - ) R Q rl + 0{r z ). 



(4.16) 
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We may now use the form of the mode functions 
to study the contribution of the u> = modes to the mode 
sum E (I2.77[) near the horizon. In other words, we wish 
to analyze 



rr, OO 

-V 

4-77- ' 

£=0 
00 

E 

1 

rfl 



4tt 

T 

-in 



(21 + 1) C MPM {r)q M {r) - —j- 

r/2 



r2! + 11 (s&)^( r (HiAu(r(/)i 



(4.17) 



where /j, B and y t are constants and x — r — r . The 
fact that this expansion should be uniform in I means 
that the constants yi must be bounded as i — » 00. The 
expansions (|4.12p lead to the following behaviour near 
the horizon: 



'"O V /l 



r(r) 

- 2 ^(r) 



■0(x). 



(4.13) 



In section HTdI we added and subtracted J2eLo flioe(r) to 
this sum to give a sum which was more readily computed 
numerically. We do not need to do this here as the sum 
is convergent. In addition, (3\qi(t) diverges as r — > tq, so 
that the sum (|4.17[) is more straightforward to analyze, 
as it does not include these additional divergences (which 
do, however, cancel). 

Using the Watson-Sommerfeld formula (|2.55p . we con- 
vert the sum (|4. 1 7[) over i into two integrals: 



The expansions (|4.12p are substituted, in turn, into 
P~TT]) . then P~TU)) and finally (|4T5j) to give the form of 



T 

An 



[25 



(4.18) 
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where 



dX 



2A 



r(r) 



I (T(r))K (T(r)) 



yll _ 



1 

-1? 



(4.19) 



dX- 



a 2tr\ 



2iX 



Io(r(r))tfo(T(r)) - — - 
rf 2 



(4.20) 



In (|4~T9|) we have I = X - 1/2, while in (|OD|) we have 
I = iX — 1/2. In both cases the functions T(r) and fi(r) 
depend on I. 

We now use the expansions (|4.13[) for the functions 
r(r) and fi(r) near the horizon, which are uniform in I. 
It turns out that Eg 7 is the more straightforward to work 
with, so we examine this first. Using (|4.13[) . the leading 
order behaviour of E 7 / near the horizon is given by 



■Jl 



3i 



dA- 



8A 



^o 2 /i(l 

25 

+O(a;ln[x/ro]) 



x/ 



2B 
ro 



(4.21) 



where -B depends on A. In this case we can use the ap- 
proximate behaviour of the modified Bessel functions for 
small x (3~l| : 

I {x) = l + 0(x 2 ); 

K (x) = -C-ln(|) + 0(x 2 \nx); (4.22) 

where C is Euler's constant. We then obtain the following 
expression for Eq 7 : 



yll _ 



1 



C + ln 
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• In 
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,2ttA 



\nB 



(4.23) 



The integral in (I4.23|l involving In B cannot be computed 
exactly because we do not have a closed-form expression 
for B for all values of £. However, this integral does not 
depend on x, and is finite because the behaviour of B for 
large A is given by ()4.16|) . Therefore we have, near the 
horizon, 



yll 



1 



In 



finite terms. 



(4.24) 



Next we consider E„ (14. 19|) . From (|4.14|) with 



we have 



2B 



dB_ 

dA 



2X 



dyi 
dX ' 



(4.25) 



so that, for large A, 



2B^- = 2X + 0(X- 2 ). 
dX 



(4.26) 



Defining Bq = B(X = 0), and using the expansions 
(I4.13p . the leading order behaviour of E 7 , can be writ- 
ten as 



E - El 



0(x ln[x/r ]), 



(4.27) 



where 

Ei = 

E 2 = 



, f AB 2B fx\ 2B fx 



dX 
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dB roVfix 
4 



dX 



X — B- 



oJi ° V 7 'o V h) ° V r o V h 



dX 



(4.28) 



To compute Ei, we first change the upper limit of the 
integral from A = oo to A = Ai > 1, which simplifies 
the analysis. The integral Si is in fact regular as Xl — > 
oo, but the individual terms in it are not. We use the 
standard result [23|, valid for all non-zero v. 



2BI Q (Bis)K (Bv)dB = B 2 [I (Bu)K (Bu) 



-hiB^K^Bv)} 



(4.29) 



which gives 



Ei 



B Jo 



2B l~x 



ro 



2B fx 



JiJ \ ro V h 



2B pF 
ro V h 



X(B) 



roVfix 



B L 



(4.30) 



where Bl = B(Xl)- For large z, we then use the expan- 
sion 1311] : 



I n {z)K n {z) 



1 

2z 



1 



4n 2 



1 



0(z~ 



(4.31) 



together with the behaviour of B for large £ l|4.16p . This 
enables us to simplify the expression (|4.30|) to obtain a 
quantity which is manifestly finite as A^ — > oo, and then 
we can take the limit Xl — > oo. The answer obtained is 



Ei = - 



2B 2 



rg/i 



ro V hJ \ r o V h 



(4.32) 
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For small x, we now expand Ei, using (|4.22|) and [3l| 



7i(i)Jfi(i) = - + 0(x 2 ln[x/r }), 



to give, for small x, 



Ei = 9 ml — J + finite terms 



(4.33) 



(4.34) 



temperature, we have shown analytically that (<f) 2 ) lcn is 
finite at both the event and cosmological horizons. This 
means that the divergences seen (Fig. [§]) in ((ft 2 ) analytic 
and (</> 2 )numeric when the scalar field is no longer massless 
do in fact cancel to give a finite ((j> 2 ) r cn- 



CONCLUSIONS 



For E2 , it is sufficient to use the expansions (|4.22[) in the 
integrand to give, for small x, 



E 2 = - 



+fmite terms 



In I - 



1 



1 

12 



dX 



X In I — 



finite terms. 



R r 2 
(4.35) 



Combining (|4.34j) and (|4. 35(1 . we find the following ex- 
pression for E^ (|4~27]) : 



T 1 
^0 



12 



+ m 2 r 2 + [ £ 



Ror, 



In 



+finite terms, 



r o, 
(4.36) 



and, combining this with (|4.24[) . gives our final expression 
for the behaviour of the mode sum E (|2.T7|) near the 
horizon: 



E = 



T 

1 

16tt 2 



Inl * 

I * 

ro 



finite terms 

finite terms, 
(4.37) 



where we have used the fact that the temperature T = 
k/271", where k = fi/2 is the surface gravity of the hori- 
zon. 

From (|4.37[) we see that the mode sum E diverges near 
a horizon unless the quantum scalar field is massless and 
conformally coupled. This behaviour can be seen in our 
numerical results in Figs. [3H3 where it is apparent that 
E is regular at a horizon in the massless, conformally 
coupled, case but otherwise divergent. 

Finally, comparing (|4. 2|) and (|4.37[) . we see that the 
logarithmic divergences cancel and E + (</> 2 ) analytic is reg- 
ular at the horizon. Since we have already shown that the 
remaining contributions to (</> 2 ) re n, namely ('Z' 2 ) analytic 
(|4.2p and (0 2 )n Umeric <|4.3[) . are regular at a horizon, wc 
therefore conclude that the total (4> 2 ) lcn , for a thermal 
state at the same temperature as the horizon, is regu- 
lar at that horizon. For the particular case of lukewarm 
black holes, the analysis above applies equally well to 
the event and cosmological horizons, which are at the 
same temperature. Therefore, for a thermal state at this 



In this paper we have studied the renormalized ex- 
pectation value (0 2 )rcn for a quantum scalar field on a 
lukewarm black hole background, where the space-time 
possesses a black hole event horizon and a cosmological 
horizon which are at the same temperature. We have 
used the method of [l( to compute (</> 2 ) ren for a thermal 
quantum state at this natural temperature. Our numer- 
ical computations have indicated that (4> 2 ) le n is regular 
on both the event and cosmological horizons, and this 
has been proved analytically as well. 

The method of Ref. [l[ splits ((/> 2 ) rcn into two parts, 
an analytic expression ((ft 2 ) analytic and a part, ((ft 2 ) numeric, 
which can only be computed numerically. One might 
hope that the analytic expression would be a good ap- 
proximation to the exact quantity, at least in some limit 
(say, for large mass). However, as observed in Ref. [l[, 
this is not the case because ((ft 2 ) analytic contains a term 
which diverges logarithmically near a horizon unless the 
quantum scalar field is massless and conformally coupled. 
We have shown that, for a thermal state at the same tem- 
perature as the horizon, this divergence cancels with a di- 
vergence in (<p 2 ) numeric near the horizon to give an overall 
finite quantity. In the literature other analytic approxi- 
mations have been developed [13, HH 1 some of which do 
not have this logarithmic divergence. 

Unlike the Schwarzschild-de Sitter black hole for 
lukewarm Reissner-Nordstrom-de Sitter black holes we 
have shown that a static quantum state can be con- 
structed which has a regular renormalized expectation 
value (</) 2 ) re n across both future and past event and cos- 
mological horizons. This is the equivalent of the Hartle- 
Hawking state [l2j for these black holes. This result is 
in accordance with the theorems of Kay and Wald [l8l |. 
who proved that no regular thermal state can exist on 
a Reissner-Nordstrom-de Sitter black hole for which the 
temperatures of the event and cosmological horizons are 
not equal. When the temperatures do match, the fact 
that the region between the event and cosmological hori- 
zons has a regular Euclidean section [13] allows one to 
construct a thermal state at the natural temperature and 
then we have shown that this state has regular ((j> 2 ) rcn on 
both the event and cosmological horizons. 

Our work in this paper has considered only (</> 2 ) ren , and 
not the RSET. The computation of the RSET, following 
[l[, mirrors that here for ((f) 2 ) TC n, but is considerably more 
complex. In particular, even though we have shown that 
(</> 2 ) ren is finite on both the event and cosmological hori- 
zons, this does not guarantee that the RSET will be finite 
there. We plan to report on this in the near future. 
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